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ABSTRACT 

We have developed a simple model to investigate the modifications of the hydrodynamics and non- 
equilibrium ionization X-ray emission in young supernova remnants due to nonlinear particle acceleration. 
In nonlinear, diffusive shock acceleration, the heating of the gas to X-ray emitting temperatures is 
strongly coupled to the acceleration of cosmic ray ions. If the acceleration is efhcient and a significant 
fraction of the shock ram energy ends up in cosmic rays, compression ratios will be higher and the shocked 
temperature lower than test-particle, Rankine-Hugoniot relations predict. We illustrate how particle 
acceleration impacts the interpretation of X-ray data using the X-ray spectra of Kepler's remnant, 
observed by ASCA and RXTE. The thermal X-ray emission provides important constraints on the 
efficiency of particle acceleration, in complement to nonthermal emission. X-ray data from Chandra and 
XMM Newton, plus radio observations, will be essential to quantify nonlinear effects. 

Subject headings: Supernova remnants — acceleration of particles — hydrodynamics — X-rays: ISM 
— ISM: individual (Kepler) — ISM: abundances 



1. INTRODUCTION 

The strong forward and reverse shocks in young 
supernova remnants (SNRs) heat the ambient interstellar 
medium (ISM) and supernova ejecta material to X-ray 
emitting temperatures. This X-ray emission contains 
information on the supernova (e.g., type, explosion energy, 
mass and composition of ejecta) and on the nature of the 
ambient medium (e.g., uniform, stellar wind). 

In addition to heating the plasma, the forward and 
reverse shocks are thought to accelerate some fraction 
of the shocked material to cosmic ray energies. A large 
body of evidence suggests that this acceleration can be 
quite efficient with as much as 50% of the ram kinetic 
energy of the ejecta being placed in relativistic particles 
and re moved from the thermal plasm a (e.g., Blandf ord fc 



Eichler 1987; Jones & Ellison 1991; Kan g & Jones 1991 



Berezhko, Ksenofontov, & Petukhov 1999), and some early 
work investigate d its o bservational impac t (He avens 1984, 
Boulares & Cox 1988, Dorfi & Bohringer 1993). However, 
despite the expected efficiency of shock acceleration, the 
modeling and interpretation of X-ray observation from 
SNRs are generally done assuming that the shocks do not 
place a significant fraction of their energy in cosmic rays. 
The recent discovery of nonthermal X-ray emission in 

shell-like SNRs li ke SN 1006 (Koyama et al 

A (Allen et al. |l99^), G347.3- 5 (S lane et al 

and RCW86 (Borkowski et al. |1988D has shown that 



1995) , Cas 



1999) 



SNR shocks do accelerate electrons to TeV energies and 
supports the case for efficient proton acceleration and the 
necessity of investigating its effects in the modeling of 
thermal X-ray emission. 

Besides providing a more self-consistent SNR model, X- 
ray models of cosmic ray modified SNRs provide a test of 
the fundamental assumption that SNRs are the primary 
sources of galactic cosmic ray ions. Since shocks put 
more energy into accelerated ions than electrons, nonlinear 
effects seen in X-ray emission will be evidence for the 



efficient shock acceleration of ions as well. 

Here, we investigate the effects of efficient particle 
acceleration on the hydrodynamics and thermal X-ray 
emission of young SNRs. We have developed a mode l 
which couples self-similar hydrodynamics (Chevalier p_983 ) 
to nonl inear diffusive shock acceleration (Berezhko & 

|1999| ; Ellison, Berezhko, 



EUison 1999; Berezhko et al. 



& Baring 2000), and X-ray emission including non- 
equilibrium ionization effects (Decourchelle & Ballet 



1994). We show that a modified self-similar approach 
is re evant in most cases and compare it to the standard 
unmodified test-particle case. 

2. NONLINEAR SHOCK ACCELERATION 

There is convincing direct and indirect evidence from 
spacecraft observations and plasma simulations that 
collisionless shocks are efficient nonlinear accelerators. 
Consequently, they self-generate magnetic turbulence 
and develop a smooth shock precursor (e.g., Drury 



1983 



1987; Jones & Ellison 1991 



1999 



Blandford & Eichler 
Ac hterbe rg, Blandford, fc R eynolds |1994| ; Terasawa et 
al 



Ellison et al. 200C). When the acceleration is 



nonlinear, the shock compression ratio increases and the 
shocked temperature decreases compared to test-particle 
values. The entire proton spectrum must be accounted 
for self-consistently with a non-thermal tail connecting the 
quasi-thermal population to the energetic one; any change 
in the relativistic population must impact the thermal 
population. In contrast, the slope and normalization of 
test-particle power laws can be changed without changing 
the temperature and density of the shock heated gas. A 
description of nonlinear effect s and the model used here is 
given in Berezhko & EUison ( |1999| ), EUison et al. ( |200G| ). 



3. MODIFIED HYDRODYNAMIC MODEL 
3.1. Model 
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The evolution of young SNRs can be de scribe d by self- 
similar solutions (Clievalier |l982 ; Nadezhin 1985 ) provided 
that the initial density profiles in the ejecta and in the 
ambient medium have power law density distributions, 
respectively p — g" r~^t^~^ (with n > 5) and p — q r^'' 
(with s < 3), where r is the radius, t the SNR age 
and g and q are constant, g can be expressed in terms 
of the ex plosio n energy and ejected mass (Decourchelle 
& Ballet 1994). The effects of cosmic ray press ure on 
SNR dynamics were first considered by Chevalier ( |1983 ) 
using a two-fluid, self-similar solution with a thermal gas 
(adiabatic index 7 = 5/3) and a relativistic gas (7 — 4/3) 
having an arbitrary const ant f raction of the total energy. 
Here, we use Chevalier's ( |1983 ) solutions for driven waves 
but with boundary conditions calculated from the shock 
acceleration model from Berezhko & Ellison (1999). 

As a function of shock velocity, Vg, ambient density, po, 
unshocked magnetic field, Bq, and ambient temperature 
To, as well as the maximum particle energy, E^max and 
injection efficiency, ?7inj,p (i-e., the fraction of protons 
with superthermal energies), the nonlinear model yields 
downstream values for the density and pressures in 
the thermal and relativistic gases, at the forward and 
reverse shocks. Alfven wave heating in the precursor 
is used, reducing the efficiency compared to adiabatic 
heating. Starting with test-particle solutions at a given 
age, we iterate until the modified solutions and nonlinear 
boundary conditions are satisfied. 

Our solutions neglect cosmic ray diffusion (assuming 
they are spatially coupled to the gas - an excellent 
approximation for all but the highest energy particles), 
and assume that both the compression ratio and the 
relativistic to total pressure ratio are independent of time 
at either shock, which is required for self-similarity to be 
valid. 

Unless explicitly stated, we take Tq = 10^ K and Bq = 5 
pG at both shocks. For the SNR parameters, we take: 
g = 3.4 10^ CCS, corresponding to Es^iMQ/M^jf/^ = 
0.53 10^^ erg, and q = po = 0.42 amu/cm^. 

3.2. Validity of the self-similar approach 

We can a posteriori determine the validity of the self- 
similar assumption, by examining the variation versus time 
of the compression ratio and relativistic to total pressure 
ratio, at the reverse shock (RS) and at the forward shock 
(FS). As seen in Figure 1, these ratios are nearly constant 
for injection parameters ?7inj,p <^ 2 10~* and negligible for 
Vin},p ^ 10~^, so that the self-similar solutions are well 
justified in these cases. 

For values of ?7inj,p ^ 2x10^^^, unmodified solutions can 
occur at very high sonic Mach numbers Mg (i.e., early ages) 
which show a rapid transition to nonlinear solutions as 
Ms drops (see Berezhko & Ellison 1999 ). The self-similar 
approach is valid if this rapid transition occurs early or 
late compared to the age of the SNR. 

3.3. Results on the dynamics 

Figure 2 shows the density and temperature as a 
function of radius for various values of ?7inj,p- As r?inj,p 
increases, the compression ratio becomes higher (up to 15 
here instead of 4), while the shocked temperature drops 



(nearly a factor of 10 between the test-particle result and 
= IQ-'^). Such trends were expected from Chevalier 
but the actual nonlinear effects lead to larger 



mi.D 



( |l98^ ), 



modifications because of the influence of particle escape: 
the maximum compression ratio is no longer limited to 7, 
but can rise, in principle, to infinite values. 




1000 



Age (year) 



Fig. 1. — As a function of remnant age for different values of 
ryjjjj p, we show (top panel) the fraction of pressure in relativistic 
particles and (lower panel) the compression ratio in the forward 
shock (dashed curves) and in the reverse shock (solid curves). 



Notably, these modifications occur for 77inj,p ^ 5x lO"'^, 
considerably lower than that deduced for the Earth's bow 
shock and fr om p lasma simulations (i.e., '7inj,p ^ 10^'^; 



1993). Table 1 gives the relative fractions of 
thermal, Eth, cosmic ray, E^, and escaped 



Ellison ct al 
kinetic, i?kin, 

energy, Ecsc, summed over the shocked ejecta and the 
shocked ambient medium. An important fraction of the 
energy can be carried off by escaping particles. Note that 
in young SNRs, the total energy in the shocked region 
increases with time as more and more ejecta is being 
shocked. 

Cosmic ray modified dynamics lead to a thinner 
interaction region between the forward and reverse shocks 
and much larger density and pressure gradients. This 
suggests that the growth of hydrodynamic instabilities in 
the shocked material near the contact discontinuity will be 
faster than otherwise. With a forward shock closer to the 
contact discontinuity, this may help these turbulent effects 
extend to the forward shock and beyond. 

Equally important, the higher densities and lower tem- 
peratures will lead to more rapid Coulomb equilibration 
than expected in a test-particle case. 



^Emax is determined from the shock radius, age, and diffusive mean free path (e.g., Ellison, Berezhko, & Baring 2000). 
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Fig. 2. — Density and temperature versus radius for different 
values of the injection efficiency. 

3.4. Velocity versus shock temperature 
Table 1. Temperature- Velocity Relation and energy 



fractions 






X 


-Ekin 




Ed- 


Ecsc 


RS 




0.17 


0.85 


0.13 


0.01 


0.00 




10-'' 


0.02 


0.80 


0.02 


0.08 


0.10 




10"^ 


0.01 


0.78 


0.01 


0.08 


0.12 


FS 




0.18 


0.48 


0.48 


0.03 


0.00 




10-^ 


0.02 


0.42 


0.04 


0.26 


0.28 




10"^ 


0.01 


0.41 


0.01 


0.24 


0.33 



We emphasize that when particle acceleration is effi- 
cient, the relation between the shock velocity and the 
mean post-shock temperature of the shocked gas is no 
longer approximated by the expression for a strong test- 
particle shock, i.e., x = [l/{^mu)] {kTs/V^) ^ 3/16. Here 
fimn is the mean particle mass and Tg is the post-shock 
temperature. Table 1 shows how x drops well below 3/16 
for large values of ?7inj_p, and indicates that shock velocities 
derived from a measured Ts may be underestimated by 
a large factor. Independent measurements of Vs and 
Ts provide strong constraints on the amount of energy 
going into relativistic particles. This was shown using 
Chandra observations of E0102. 2-7219 SNR, by Hughes, 
Rakowski, & Decourchelle (2000), who concluded that the 
low observed post-shock electron temperature combined 
with the high observed shock velocity could only be 
understood if a non-negligible fraction of the shock ram 
energy went into cosmic rays. 

4. INTERPRETATION OF X-RAY THERMAL SPECTRA 

The increase in shocked density and decrease in shocked 
temperature from nonlinear effects, strongly impact the 
thermal X-ray emission. Compared to test-particle 



models, we predict that regions in young SNRs undergoing 
efficient acceleration will have a lower temperature, more 
intense thermal X-ray emission (due to the larger emission 
integral J dV and lower temperature), and a faster 
ionization timescale. In older remnants, the temperature 
should drop below X-ray emitting values earlier. 

We have coupled our modified hydrodynamics model 
to the calculation of the non-equilibrium ionization and 
thermal X-ray emission. The ionization and rccombii iation 
rates have been taken fro m Ar naud & Rothenflug (1985) 
and Arnaud & Raymond ( |1992| ). The X -ray emi ssion code 
is from Mewe, Gronenschild, & V an de n Oord ( 1985| ) and 
Mewe, Lemen, & Van den Oord ( 1986| ). The abundances 
in th e interstellar medium are taken to be solar (Meyer 
1985 ). We assume full equipartition between the electrons 
and ions (Tg = Ti) and neglect any contribution from the 
superthermal electrons to the ionization or line excitation. 
For efficient shock acceleration, this may be important and 
we are investigating these effects (Porquet et al. 2000). 
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Fig. 3. — Best-fit model to tfie spectra of Kepler's SNR. The 
contribution from the shocked ejecta (RS) is shown with dashed 
lines, the shocked ambient matter (FS) is shown with dotted lines, 
and the sum is shown with solid lines. 

As an indication of how the interpretation of X-ray 
data depends on acceleration efficiency, we model the X- 
ray observations of Kepler's SNR obtained over a broad 
energy range with ASCA (Kinugasa & Tsunemi 1999) 
and RXTE (Decourchelle & Petre 1999; Decourchelle et 
al. 2000). For several values of r?inj,p, we adjust the 
SNR parameters to fit the observed spectrum, which is 
dominated by the emission from the shocked ejecta as 
shown in Figure 3. All values of ?7inj,p yield similar values 
for the reduced chi-squared ('^ 4), the ionization timescale, 
the reverse shock temperature ~ 3 keV, the hydrogen 
column density ^ 4 10^^ cm^^, and oversolar abundances 
of heavy elements in the shocked ejecta (O, Ne, Mg, Si, 
S, Ar, Ca, and Fe). The derived values of forward and 
reverse shock compression ratios (rps and trs), FS and 
RS speeds {V-ps and Vrs), ambient density, distance, and 
the energy-mass relation are, however, very different as 
shown in Table 2. 

To get a temperature high enough to produce the iron 
K-line and the high energy continuum, all nonlinear fits 
require larger shock velocities and lower ambient densities. 
These imply more kinetic energy in the ejecta and a larger 
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Table 2. Best-fit derived quantities for different values of the injection parameter at both shocks. 



»?inj,p FS 


J^inj.p RS 


kTps 
keV 


keV 


rFS 


^■RS 


V-ps 
km/s 


^^RS 
km/s 


po 

amu/cm'^ 


D 
kpc 


iJsN(A/0/Mej)^/» 

10^^ erg 








17 


3 


4 
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4000 


1700 


0.67 


5 


0.3 


2 10"* 
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14 


13 


12200 


5800 


0.23 
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10"^ 
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11 
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17 


16 


15800 


7600 


0.18 


20 


18 


10"^ 


10"^ 


11 
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18 


17 


18200 


8800 


0.17 


23 


30 


10"^ 





1.5 
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13 


4 


4570 


2110 


0.57 
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0.5 



inferred distance, D, than the test-particle case. To obtain 
a distance consiste nt wi th observations (4.8 ±1.4 kpc, 
Reynoso and Goss 1999 ), a low efficiency at the reverse 
shock is required. This can be done in three ways: with 
Bq ^ 25/iG at the reverse shock large enough to dampen 
acceleration, with Bq ^ 0.1 and yyinj.p rs ^ lO^'' to 
allow high Mach number, unmodified solutions, or with 
no reverse shock acceleration (jyinj.p rs — 0)- A low 
magnetic field is expected from the expansion of the ejecta, 
but magnetic field amplification could increase it. In all 
cases, a good fit is obtained for efficient acceleration at the 
forward shock and test-particle conditions at the reverse 
shock. Figure 3 and Table 2 (last row) show the case with 

'yinj.p RS = 0- 

Spatially resolved spectroscopy is required to confirm 
both a low efficiency of the acceleration at the reverse 
shock and our prediction of efficient acceleration at the 
forward shock, which is expected for standard values of 
the upstream density and magnetic field. The arbitrary 
power law component required to fit the broad-band X-ray 
continuum (Decourchelle & Petre 1999 ) is likely related to 
the forward shock. In future work, this component will 
be replaced by including self-consistently the nonthermal 
continuum, and using additional constraints on the 
acceleration parameters from the radio data. 

5. CONCLUSION 

Nonlinear effects from efficient particle acceleration 
strongly influence the X-ray emission in SNRs, and may 
serve as a tracer of cosmic ray acceleration. Here, we 
have focused on the consequences of shock acceleration 



on the thermal X-ray emission and have presented the 
first nonlinear, non-equilibrium ionization modeling of 
the X-ray emission from a specific SNR. While our self- 
similar approach has limitations which must be carefully 
considered, we believe it contains the essential physics and 
allows the exploration of a large range of parameters. We 
find that typical source parameters of young SNRs should 
result in nonlinear acceleration at the forward shock, but 
yield near test-particle conditions at the reverse shock if 
the magnetic field is as low as expected from the ejecta 
expansion. 

When cosmic ray production is efficient, shock heating 
is linked to particle acceleration and thus to the broad- 
band photon emission from energetic electrons and ions. 
The high energy electrons produce both radio and X-ray 
synchrotron and this will be included in future work. 

In addition to a possible nonthermal synchrotron 
component, regions undergoing strong acceleration should 
exhibit stronger and lower temperature thermal emission, 
while regions with little acceleration should be fainter and 
at a higher temperature. 

Also, other complications (e.g., the vari ation of the 
ambient magnetic field orientation (Reynolds |l998 ), shock 
interactions with clouds and dense knots, etc.) influence 
the full modeling of SNRs. The spatially resolved spectra 
from the Chandra and XMM-Newton satellites should 
greatly aid the investigation of these and other effects. 

We thank E. Berezhko, L. Drury, J. Hughes, J. P. Meyer 
and E. Parizot for fruitful discussions as well as the referee 
for helpful comments. 
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